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1.  Introduction 


There  now  exists  substantial  observational  evidence  for  strong  heating 
of  the  temperature-minimum  )  region  in  solar  flares  (Machado  et  al. . 

1978;  Cook,  1979;  Cook  and  Brueckner,  1979;  Canfield  et  al.,  1980).  This 

observed  temperature  rise  amounts  to  some  100-200  K,  corresponding  to  a 

-3  -1 

specific  energy  input  of  some  10-20  erg  cm  s  at  these  levels  (Machado 
et  al.,  1978).  It  appears  to  be  difficult  to  reconcile  this  amount  of  heat¬ 
ing  with  a  primary  flare  energy  release  in  the  corona  and  subsequent  transport 
of  this  energy  through  the  atmosphere  by  either  particles  or  radiation 
(Machado  et  al. ,  1978;  Emslie  and  Machado,  1979).  Emslie  and  Machado  (1979), 

among  others  (e.g.,  Sturrock,  1980),  have  therefore  proposed  that  the  T  . 

min 

regioa  is  heated  by  joule  dissipation  of  currents.  They  do  not  consider 

the  generation  of  these  currents,  but  they  do  point  out  that  such  currents 

must  have  a  very  fine  spatial  structure  (£  106  cm)  since  otherwise  the 

magnetic  field  associated  with  them  would  be  unacceptably  large. 

According  to  Machado  et  al.  (1978).  the  total  energy  dissipated  at  T  . 

— - -  min 

levels  during  flares  is  small  compared  to  the  total  flare  energy  budget.  Thus, 

in  explaining  T  heating,  the  total  energy  release  requirement  is  not 
min 

difficult  to  meet.  The  problem  instead  is  to  find  a  mechanism  for  transport¬ 
ing  this  energy  efficiently  from  the  primary  release  site  to  the  T  ,  level 

min 

of  the  atmosphere  (or  to  explain  how  part  of  the  energy  can  be  released  at  the 
Tmin  leVGl)* 

Emslie  and  Machado  (1979)  note  that,  because  of  the  relatively  low  tempera¬ 
tures  and  low  fractional  ionization  level  of  the  plasma  at  T,  levels,  the 

min 

plasma  resistivity  r\  there  is  much  larger  than  at  any  other  point  in  the 
solar  atmosphere  (see  Figure  1).  Hence  the  energy  carried  by  electric  currents 
generated  in  the  corona  may  reach  the  T^^  layers  of  the  atmosphere,  essen¬ 
tially  unattenuated,  and  then  be  efficiently  damped  due  to  the  large  resistivity. 
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The  transverse  scale  of  such  coronal  current  patterns  generated  during 
flares  may  well  be  small  enough  (see  Section  3)  to  satisfy  the  constraints 
imposed  by  Emslie  and  Machado  (1979). 

In  this  paper  we  therefore  Investigate  the  possibility  that  electric 
currents  associated  with  Alfven-wave  perturbations,  generated  by  the  primary 
flare  energy  release  in  the  corona,  are  responsible  for  heating  the  tempera¬ 
ture-minimum  region  during  flares.  In  Section  2  we  discuss  the  propagation 
and  damping  of  these  waves.  In  Section  3  we  investigate  the  heating  rates 
effected  by  this  mechanism  for  a  range  of  model  parameters  and  so  place 
constraints  on  these  parameters  if  the  mechanism  is  to  be  effective  at  T 
heating.  In  Section  4  we  discuss  the  results  obtained. 

2.  Propagation  and  Damping  of  Alfven  Waves  in  the  Solar  Atmosphere 

We  wish  to  analyze  the  transmission  of  energy  along  magnetic  field  lines 
from  coronal  heights  down  to  the  temperature-minimum  region.  In  the  MHD 
approximation,  there  are  three  modes  of  propagation:  the  slow  and  the  fast 
mangetoacoustic  modes,  and  the  Alfven  mode.  In  a  complex  situation  involving 
non-uniform  plasma  and  magnetic  field,  the  group  velocity  of  the  magneto¬ 
acoustic  modes  must  be  expected  in  the  course  of  propagation  to  depart  from 
the  direction  of  the  magnetic  field,  even  if  the  group  velocity  is  initially 
directed  along  the  field  (Weinberg,  1962).  On  the  other  hand,  the  group 
velocity  of  Alfven  waves  is  always  parallel  to  the  magnetic  field.  Energy  in 
the  Alfv&i  mode  will  therefore  remain  in  that  mode  as  long  as  the  wavelength 
is  small  compared  with  the  length  scale  for  variations  of  the  field  and  plasma. 

The  scale  height  for  variation  of  the  plasma  properties  varies  over  a 
wide  range  in  the  sun's  atmosphere.  The  wavelength  of  a  propagating  Alfven 
wave  also  varies  with  height  since  the  Alfvln  speed  varies.  We  find  that,  for 


waves  of  interest,  the  wavelength  is  much  less  than  the  scale  height  for 
variation  of  density  (for  instance)  throughout  the  corona  and  in  the  upper  part  of 
the  transition  region  (see  Section  3).  We  begin  by  considering  this  region 
and  postpone  until  later  in  this  section  discussion  of  the  lower  transition 
region  in  which  the  density  varies  abruptly. 

Since  we  are  concerned  only  with  the  propagation  of  energy  along  magnetic 
field  lines,  we  shall  ignore  the  magnetoacoustic  modes  and  consider  only  the 
Alfven  mode.  We  begin  by  considering  a  simple  model  in  which  the  magnetic 
field  is  uniform  and  vertical,  but  we  shall  allow  for  non-zero  resisitivity 
and  for  variation  of  density  with  height.  We  shall  subsequently  consider  the 
additional  effect  of  slow  variation  of  the  magnetic  field  strength  with  height. 

The  equation  of  motion  of  the  plasma  is 

PdF="^+^  *  (2-1) 

where  we  use  modified  Gaussian  units  (Jackson,  1962)  for  electromagnetic 

quantities,  and  we  have  ignored  the  gravitational  force  since  we  shall  be 

concerned  with  wave  frequencies  much  higher  than  the  Brunt-Vaisala  frequency 

(Unno  et  al. ,  1979) .  The  unperturbed  state  is  characterized  by  the  density 

p  (z)  ,  the  pressure  p  (z)  and  magnetic  field  strength  B  (z)  «  (0,0, B  )  . 
o  o  o 

To  first  order,  the  perturbation  in  the  equilibrium  state  of  the  plasma  and 
field  satisfies 

d6y 

P0  - - £5p  +  5^o  .  (2.2) 

On  using  the  Maxwell  equations 

7  (2‘3) 
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and 


(2.4) 


and  che  Ohm's  law  relationship 

^4  *  i +  i  ^  * 


where  the  resistivity  n  is  taken  to  be  a  scalar,  we  obtain 


dl 

3t 


+  4?  ^  -4t^W  ‘ 


(2.5) 


(2.6) 


On  combining  equations  (2.2),  (2.4)  and  the  appropriate  perturbed  form  of 
(2.6),  we  obtain  the  wave  equation 


at* 


{["  ^  (*5p)  +  4 sr  w  x  «o]x 


(2.7) 


Since  the  Alfven  mode  is  transverse,  we  may  without  loss  of  generality  consider 
that  *  (0,  6By ,  0)  .  Further,  since  the  AlfvSa  mode  involves  divergence- 
free  motion,  the  pressure  gradient  plays  no  role.  (Our  assumption  that  the 
wavelength  is  small  compared  with  the  scale  height  in  fact  makes  it  possible 
to  neglect  this  term  not  only  for  the  Alfv£n  mode  but  also  for  the  fast  and 
slow  magnetosonic  modes.)  Hence  we  may  write 


6B 

y 


b(z)e 


i(kxx  -  tot) 


(2.8) 
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Equation  (2.7)  then  reduces  to 


H  +  J)  b  • 

where  V  is  the  Alfven  speed  given  by 


(2.9) 


V2 


4np 


(2.10) 


If  (hence  V)  is  independent  of  z  and  if  n  *  0  ,  then  b 

varies  as  e^z2,  where  k  satisfies  the  dispersion  relation  for  Alfven 

z 

waves : 

0)2  -  V2  k2  .  (2.11) 

z 


Since,  for  the  problem  of  interest,  n  #  0  >  the  Alfven  waves  will  be  damped. 
If  V  varies  slowly  with  z  ,  then  an  Alfven  wave  will  retain  its  identity 
but  its  amplitude  will  vary  in  the  course  of  propagation.  We  may  consider 
both  of  these  effects  at  the  same  time  by  expressing  b(z)  in  the  form 


b(z) 


A(z)  exp 


z 


k 

z 


(z ’ )dz 


(2.12) 


We  find  that,  to  lowest  order  in  the  effects  of  n  and  the  variation  of  V  , 
the  expression  given  by  (2.12)  indeed  satisfies  the  wave  equation  (2.9)  if 
the  amplitude  A(z)  varies  as  follows: 


A(z)  =  A(0) 


(2.13) 


In  deriving  this  expression,  we  have  assumed  that  the  wave  is  propagating  in 
the  z-direction.  The  "damping  length"  is  expressible  as 


6 


(2.14) 


where  and  L„  ,  the  "transverse"  and  "parallel"  damping  lengths,  are 

given  by  the  expressions 


Ll 


4ttV  _  4uV  l2 
cnk^  cn 


(2.15) 


and 


t  4irV3 
11  *  cnw2 


(2.16) 


and  we  write  1  for  (k  )  1  . 

x 

If  the  energy  flux  being  carried  by  an  Alfven  wave  is  represented  by 
S  ,  and  if  the  actual  perturbation  in  the  magnetic  field  due  to  the  wave  is 
given  by  the  real  part  of  expression  (2.12),  then 


S 


(2.17) 


Hence,  from  Equation  (2.13), 


S(z) 


S(0)  exp 


'  r*  dz± _ ' 

V  OPT 

u  n  u  J 


(2.18) 


This  equation  will  be  used  to  estimate  the  damping  of  Alfven  waves  in  the 
neighborhood  of  the  temperature-minimum  region. 

As  an  Alfven  wave  of  given  frequency  propagates  down  through  the  corona 
into  the  transition  region,  it  encounters  a  density  gradient  which  is 
initially  small  but  becomes  larger  and  larger.  At  some  location,  the  assump¬ 
tion  that  the  wave  length  is  small  compared  with  the  characteristic  length 
scale  of  variation  of  density  becomes  invalid.  Below  that  point,  the  density 
gradient  becomes  very  steep,  so  that  there  is  in  fact  a  sharp  change  in 
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density  in  less  than  one  wavelength.  In  this  situation,  part  of  the  wave 
energy  will  be  transmitted  and  part  will  be  reflected. 

The  component  of  the  wave  vector  parallel  to  the  z-axis  is  given  by 


k 

z 


-1  -1  1/2 
<u  V  =•  oiB  (4irm  n) 

P 


(2.19) 


where  n  is  the  number  density  both  of  electrons  and  protons.  To  a  fair 

approximation,  the  transition  region  may  be  regarded  as  a  region  in  which 

~2  —1 

the  heat  flux  F  (erg  cm  s  )  is  constant  (Athay,  1976).  Then,  using 
the  Spitzer  (1962)  conductivity,  the  characteristic  length  scale  for  the 
temperature  change  L^,  ,  defined  by 

Lt  -  T  (dT/dz)"1  ,  (2.20) 

is  given  by 

Lt  -  otT 7/2  F"1  ,  (2.21) 

where  a  %  10-6  erg  cm  1  s  *K  7//2  .  The  gas  pressure,  given  by 

P  =  2nkT  ,  (2.22) 

is  approximately  constant  in  the  transition  region  so  that  the  density 
scale  height  is  equal  to  the  temperature  scale  height  L^,  .  Hence  the 
critical  condition  that  the  density  gradient  is  comparable  with  the  Alfvgn 
wave  number  may  be  expressed  as 


k  L_  »  1  .  (2.23) 

z  r 

On  using  Equations  (2.19),  (2.21)  and  (2.22),  we  find  that  this  condition 
occurs  at  the  "critical"  layer  where  the  density  and  temperature  are  n£ 
and  T£  ,  given  by 


8 


and 


n  «  1012*4  F-1/3  P7/,S  B-1  co1/3 

c 


T  *  103 ’2  F1/3  p-*1 /6  B1/3  co  1^3 
c 


(2.24) 


(2.25) 


As  far  as  transmission  and  reflection  are  concerned,  the  Important  quantity 
is  the  jump  in  density  in  going  from  this  critical  region  into  the  chromosphere, 
which  we  write  as 


p  .  n  .  T 

ch  ^  ch  m  __c 

P  n  T  , 

Kc  c  ch 


(2.26) 


We  see  from  Equation  (2.25)  that  this  may  be  expressed  as 


0  .  103-2  T-1  F1/3  P-1/e  B1/3  of1/3  . 
ch 


(2.27) 


We  write  the  wave  amplitude  on  the  coronal  side  of  the  transition  region  as 

i[k_(z-z  )  -  cot]  i[-  k  (z-z  )  -  cot] 

b_  (z,t)  ■  Ae  +  Rg  Ae  _  (2.28) 


and  the  wave  amplitude  on  the  chromospheric  (transmitted)  side  as 


b+  (z,t)  »  Ae 


i[k+(z-Zj)  -  cot] 


(2.29) 


where  the  jump  is  assumed  to  occur  at  z  »  Zj  .  On  noting  that  b  is  con¬ 
tinuous  across  the  jump  and,  from  Equation  (2.9),  that 


db  db, 

V2  - -  a  V2  - 

-  dz  +  dz 


(2.30) 


we  find  that  the  reflection  and  transmission  coefficients  are  given  by 
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„  9l/2  ~  1  T  _  2  9l/z 

0l/2  +  i  ’  B  01(,2+1 


(2.31) 


On  noting  that  the  energy  flux  is  proportional  to  b2V  ,  we  see  that  the 
coefficients  for  the  reflection  and  transmission  of  energy  are  given  by 


.  (§j^  -  I)’  .  46|/2 

h  (3i/2  +  i)2  ’  E  (e1/2  +  i) 2 


(2.32) 


If  energy  is  suddenly  released  at  the  top  of  a  loop,  then  a  fraction  T£ 
of  that  energy  is  transmitted  to  the  chromosphere  and  photosphere  after  a  time 
tj  ,  where  t^  is  the  time  it  takes  for  an  Alfven  wave  to  propagate  from  the 
top  of  the  loop  to  the  base.  After  an  additional  time  2tj ,  an  additional 
fraction  R^T^,  is  transmitted,  etc.  (We  have  assumed  that  the  wave  is  not 
significantly  damped  in  the  corona;  this  will  be  verified  a  posteriori  in  Section 
3.)  Hence  we  see  that  the  characteristic  decay  time  x^  of  the  energy  flow  to 
the  lower  atmosphere  is  given  by 


x 


D 


(2.33) 


For  example,  consider  the  range  of  values  0  =  10  -  100  .  For  this  range, 
we  find  that  R^  *  .27  -  .67  ,  so  that  Xp/tj  *  1*5  -  5.0  .  In  a  symmetrical 
configuration,  the  energy  will  escape  equally  through  both  feet  of  the  loop. 

We  have  so  far  considered  that  the  magnetic  field  is  uniform  and  vertical. 
In  a  real  solar  flux  tube,  neither  of  these  conditions  will  be  satisfied.  Wein¬ 
berg  (1962)  has  discussed  very  generally  the  effect  upon  wave  propagation  of 
variations  in  the  parameters  entering  the  wave  equation.  His  analysis  confirms 
that  for  Alfven  waves,  even  in  a  general  magnetic-field  configuration,  energy 
propagates  along  magnetic  field  lines  and  that  the  wave  amplitude  varies  accord¬ 
ing  to  Equation  (2.13). 
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However,  Weinberg's  analysis  also  sheds  light  on  the  variation  of  the 
wave  vector  during  propagation.  If  we  consider  propagation  along  a  field  line 
which  is  curved  but  lies  in  a  plane,  we  find  that  the  components  of  the  wave 
vector  along  the  field  line  (k,,)  and  normal  to  the  plane  (k^)  are  un¬ 
affected  by  the  curvature,  but  the  rate  of  change  of  the  component  k^ 
normal  to  the  field  line  acquires  the  following  contribution 


dkx 

—  =  -  2k  k,,  ,  (2.34) 


where  k(s)  is  the  curvature  of  the  field  line.  In  order  to  assess  the  influ¬ 
ence  of  this  effect,  we  may  consider  the  case  that  a  wave  begins  at  the  top  of 
a  loop  (s  =  0)  with  kA  =  0  .  Then  at  the  base  of  the  loop,  where  s  =  L  , 


ML) 

k„(L) 


L  ds  ic(s) 
[V(s)/V(L)] 


(2.35) 


As  a  result  of  the  rapid  increase  in  density  as  s  approaches  L  ,  the  dominant 
contribution  to  the  integral  comes  from  the  last  (the  lowest)  density  scale 
height.  If,  as  is  likely,  the  angular  deflection  ,/K(s)ds  of  a  magnetic 
field  line  is  small  in  the  lowest  density  scale  height,  we  may  neglect  the 
change  in  the  wave  vector  due  to  curvature  of  the  field  lines. 

The  wave  vector  will  also  be  modified  during  propagation  by  the  convergence 
or  divergence  of  the  field  lines.  If  the  convergence  is  "one-dimensional"  in 
the  sense  that  there  is  convergence  or  divergence  in  one  of  the  transverse 
directions  (x  ,  say)  but  not  in  the  other  direction  (y)  ,  then  kx  scales 
with  the  magnetic  field  strength  but  k^  is  unaffected.  If,  as  another  example, 

the  field  configuration  has  cylindrical  symmetry,  then  both  of  the  transverse 

\  I  o 

components  of  the  wave  vector  scale  as  B  .  We  shall,  for  definiteness, 
consider  the  latter  case  in  the  examples  to  be  discussed  in  the  next  section. 
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Table  1 


Figure  1 
Figure  2 


3.  Model  Parameters 

The  variation  of  b(z)  with  z  ,  as  given  by  Equations  (2.13),  (2.28) 

and  (2.29),  depends  not  only  on  the  values  of  the  wave  parameters  to  and 

£q(=  l  at  the  top  of  the  loop),  but  also  on  the  features  of  the  flare  loop 

itself.  We  have  investigated  four  possible  loop  models,  shown  in  Table  I« 

These  models  differ  only  in  their  magnetic  field  strengths  and  geometries, 

with  the  coronal  field  strengths  in  the  range  100-300  G  and  the  field  strengths 

near  the  footpoints  of  the  loop  in  the  range  100-900  G.  We  consider  a  narrow 

flux  tube  with  "cylindrical”  symmetry.  For  the  purposes  of  calculation,  a 

convenient  form  for  the  variation  of  B  with  z  is  the  model  considered  by 

o 

Antiochos  and  Sturrock  (1976), 

B°(Z)  1[z  -l/2\ 

FToT  *  sec  \l  arc  cos  r  /  ’  (3a) 

o 

(although  this  model  was  derived  for  a  field  configuration  of  somewhat  differ¬ 
ent  symmetry) .  The  area  factor  ("compression  ratio")  T  for  each  model  con¬ 
sidered  is  tabulated  in  Table  I.  Below  z  =  L  ,  the  background  field  strength 
was  taken  to  be  uniform  and  equal  to  Bq(L)  .  The  chromospheric  variations  of 
temperature,  density,  ionization  level  and  resistivity  with  height  were 

taken  from  model  'FI'  of  Machado  et  al.  (1980)  and  are  shown  in  Figure  1, 

—2 

plotted  against  the  column  number  density  N  (cm  )  as  independent  variable. 

The  resistivity  was  assumed  classical  (Spitzer,  1962),  with  allowance  made  for 
the  effect  of  ion-neutral  collisions  near  the  temperature  minimum  (Emslie  and 
Machado,  1979).  The  possibility  of  anomalous  resistivity  in  the  region  of 
flare  energy  release  will  be  considered  below. 

We  solved  the  relevant  equations  for  each  atmospheric  model  for  a  range 
of  both  u("  a>/2ir)  and  ;  the  resulting  behavior  of  the  quantity  Q  (local 
energy  deposition  rate  per  unit  volume)  is  shown  in  Figure  2.  The  results  in 
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Che  figure  correspond  Co  an  inicial  perturbation  amplitude  b(0)  -  1G  ;  for 
ocher  values  of  b(0)  ,  Che  energy  deposicion  race  scales  as  [b(0)]2  .  The 

shaded  areas  in  Che  figure  represent  Che  approximate  location  of  the  strong 
T  ,  temperature  enhancements  reported  by  Machado  et  al.  (1978);  the  maximum 
of  Q(N)  must  lie  in  this  region  if  Che  heating  mechanism  proposed  in  the 
present  paper  is  to  be  effective.  For  l  values  much  larger  than  106  cm  , 
the  L,,  term  in  Equation  (2.14)  completely  dominates  [cf.  Equations  (2.15)  and 
(2.16)], and  so  the  results  are  similar  to  those  for  »  10s  cm  .  Similarly, 
for  v  <<  1  Hz  ,  the  term  dominates,  so  that  results  for  such  low  frequen¬ 

cies  need  not  be  shown  either. 

It  is  immediately  obvious  from  the  figure  that  the  proposed  mechanism  is 
effective  at  heating  T^n  layers  only  if  the  parameters  v  and  l  fall 
within  the  ranges 

V  t  10  Hz  (3.2) 

and 

*0  Jfc  10s  cm  .  (3.3) 

Further,  only  for  models  with  high  photospheric  magnetic  field  strengths  [result¬ 
ing  in  higher  Alfven  velocities  and  so  lower  damping  rates  in  the  chromosphere 
-  see  Equations  (2.15)  and  (2.16)]  Is  a  satisfactory  situation  found  (e.g. 
models  c,  d  In  the  case  l  ■  105  cm  ,  v  -  1  Hz)  .  This  result  is  encouraging, 
since  the  results  of  Tarbell  et  al.  (1979)  Indicate  that  such  high  values  of 
the  magnetic  field  strength  at  the  photosphere  are  indeed  appropriate. 

The  required  values  of  v  and  l  may  be  compared  with  independent  esti¬ 
mates  of  their  likely  values.  Recent  analyses  of  the  time  structure  of  flare- 
associated  hard  X-ray  bursts  (Dennis  et  al.,  1981)  reveal  structures  with  time- 
scales  of  the  order  of  a  few  tenths  of  a  second,  with  no  significant  structure 
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on  cime-scales  smaller  than  this.  Thus,  based  on  this  data,  we  would  expect  the 
highest  frequency  of  any  flare-associated  disturbance  to  be  of  order  10  Hz  , 
in  good  agreement  with  the  values  required  (Figure  2).  Further,  on  the  basis 
that  magnetic  reconnection  in  an  unstable  arch  is  the  source  of  the  primary  energy 
release,  we  should  expect  %  to  be  comparable  to  the  width  of  a  magnetic 
"island"  (Spicer,  1976,  1977).  From  Spicer's  (1976)  Equations  (4.35)  and 

(4.54)  (see  also  Furth  etal . ,  1963),  we  find  that 


where  K  is  the  wavenumber  of  the  perturbation  inducing  the  magnetic  reconnec¬ 
tion.  Inserting  numerical  values  and  noting  that  K  %  L  1  ,  we  obtain,  using 

the  coronal  model  given  in  Table  I, 

li  %  1°S‘S(tsF")1/5  %  Io3'3  cm  (3,5) 

(see  also  Van  Hoven,  1976).  In  deriving  this  value  we  have  used  the  classical 
plasma  resistivity  (Spitzer,  1962).  The  presence  of  turbulence  in  the  reconnec¬ 
tion  region  will  increase  h  and  so  .  For  instance,  is  increased  to 

10s  cm  if  the  turbulent  resistivity  exceeds  the  Spitzer  resistivity  by  a 
factor  of  about  2x10**  .  Hence  the  required  value  ^  105  Cm  is  not  incom¬ 
patible  with  the  "island  thickness"  likely  to  arise  during  reconnection. 

As  mentioned  in  Section  1,  the  observations  and  empirical  model-fitting  of 
Machado  et  al.  (1978)  indicate  flare-associated  temperature  rises  of  some  100- 

200  K  at  T  levels,  corresponding  [cf.  their  Equations  (18)  and  (19)]  to 
min 

-3 

a  specific  energy  depositicn  rate  of  some  10-20  erg  cm  s  ,  depending  on 
whether  the  observed  structure  is  "average"  or  "bright".  From  the  results  of 
Figure  2  (which,  we  recall,  are  normalized  to  unit  initial  perturbation  ampli- 
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tude  b(0))  ,  we  therefore  deduce  that  the  required  value  of  the  magnetic- 
field  fluctuation  generated  by  reconnection  should  have  the  value 

b(0)  *  5G  (3.6) 

in  order  for  Q  to  agree  with  these  observationally  deduced  values.  This 
value  of  b(0)  also  corresponds  to  a  total  (i.e.  depth-integrated)  energy 
deposition  rate  of  some  5xl07  erg  cm  2s  1  in  the  layers  (see  dashed 

curves  in  Figure  2);  this  is  also  in  good  agreement  with  Machado  et  al. 's 
(1978)  estimates. 

Such  a  value  of  b(0)  results  in  a  value  of  b  of  %  15  G  just  below 
the  transition  zone  [see  Equation  (2.31)].  The  value  of  b  further  increases 
to  %  40  G  at  a  column  depth  N  %  3xl022  cm  2  ,  because  the  increase  of  p  and 
the  corresponding  decrease  of  V  [Equation  (2.10)],  leads  to  an  increase  of  b 
(Equation  (2.13)].  For  larger  N  ,  b  falls  off  due  to  resistive  damping  of 
the  waves  which  overcomes  any  increase  in  b  due  to  the  effect  of  increasing 
density.  These  values  of  b  are  much  smaller  than  the  ambient  field  strength 
Bq  in  the  models  of  interest,  as  is  necessary  to  justify  the  linear  analysis  of 
Section  2.  Finally,  we  note  that,  for  the  values  of  v  and  t,Q  emerging  from 
our  study,  the  wavelength  in  the  corona  is  much  less  than  the  temperature  length 
scale  in  that  region.  Also,  both  Ln  and  [Equations  (2.15)  and  (2.16)]  are, 
for  classical  resistivity,  much  larger  than  L  ,  so  that  the  wave  damping  in  the 
corona  (z  <  Zj)  is  negligible.  In  fact,  even  if  the  resistivity  in  the  region 
of  magnetic  reconnection  exceeds  the  Spitzer  (1962)  value  by  the  factor  of  %  10** 
required  to  reconcile  &o  with  the  thickness  of  the  reconnection  region  (see 
above),  it  is  still  easily  demonstrated  that  both  L„  and  are  large  compared 
to  the  size  of  such  a  region.  We  thus  see  that  the  modeling  assumptions  of 
X  «  L*,  and  L_  »  L  made  in  Section  2  are  justified. 
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4.  Discussion  and  Conclusions 

As  pointed  out  above  (Section  1),  the  total  energy  liberated  at  Tmin 
levels  during  solar  flares  is  a  relatively  small  factor  in  the  total  energy 
budget  of  a  flare,  and  therefore  the  constraint  that  the  observed  strong  heat¬ 
ing  at  T  .  levels  places  on  a  candidate  heating  mechanism  does  not  hinge 
min 

upon  the  total  energy  supply  but  simply  upon  the  existence  of  a  mechanism  to 
transport  energy  throughout  the  solar  atmosphere  from  the  corona  until  it  is 
efficiently  damped  at  levels.  The  mechanism  proposed  in  the  present 

paper  permits  virtually  undamped  propagation  of  energy  through  the  corona  and 
upper  chromosphere,  coupled  with  an  efficient  damping  mechanism  at 
levels. 

Ue  conclude  from  the  results  of  the  previous  section  that  Alfven  waves 

of  moderate  amplitude  can,  through  resistive  dissipation  of  their  associated 

alternating  currents,  account  for  the  observed  temperature  enhancements  at 

T  .  levels  during  solar  flares.  However,  it  is  interesting  to  note  that 
min 

the  propagation  of  energy  from  the  primary  release  site  in  the  corona  to  the 

T  .  level  is  relatively  slow,  with  an  associated  timescale 
min 

T  *  J Hz)  '  (4,1) 


For  the  coronal  and  chromospheric  models  adopted  in  this  paper  (Table  I),  we 
can  readily  compute  this  integral  numerically;  we  find  that  (with  T  measured 
in  seconds) 


T 


-  104 

B  . 
o,ph 


(4.2) 


where  B0>ph (gauss)  is  the  strength  of  the  photospheric  magnetic  field  [the  exact 
field  geometry  in  the  corona  is  unlmportane,  since  the  disturbance  travel  time 
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in  the  corona  represents  only  a  small  contribution  to  t;  see  Equation  (2.33)]. 

This  is  much  larger  chan  the  time-scale  quoted  by  Emslie  et  al.  (1978)  for  an 

MHD  disturbance  to  propagate  in  the  upper  chromosphere,  because  of  the  much 

larger  densities  (and  correspondingly  low  Alfven  velocities)  appropriate  to 

the  T  .  region.  Nevertheless,  for  the  plausible  range  of  values  B  , 
min  o , pn 

•i  300  G  ,  the  timescale  is  much  smaller  than  the  10  minute  delay  time  between 
the  impulsive  phase  of  the  flare  and  the  T^  enhancements  observed  by 
Machado  et  al.  (1978).  This  therefore  implies  that,  for  the  proposed  heating 
model  to  be  valid,  magnetic  disturbances  must  be  generated  in  the  corona  (or 
elsewhere)  long  after  the  impulsive  phase  of  the  flare.  It  is  possible  that 
these  disturbances  can  reveal  themselves  in  other  observational  signatures, 
such  as  radio  emission;  however,  it  is  unlikely  that  their  associated  currents 
will  produce  observationally  detectable  effects  in  locations  other  than  the 
Tmin  re8i°n»  because  of  the  very  low  resistivity  in  all  other  parts  of  the 
atmosphere. 

Finally,  we  note  that  Equation  (4.2)  represents  a  strong  observational 

test  of  the  present  model.  If  Equation  (4.2)  is  valid,  no  enhancement  at  the 

T  .  level  should  be  detected  until  a  time  measurably  later  than  the  impulsive 
min 

phase  of  the  flare  —  as  indicated,  for  example,  by  the  hard  X-ray  burst.  It 
appears  that  data  do  not  yet  exist  against  which  to  check  this  consequence 
of  the  model. 
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Figure  Captions 


Figure  1.  Chromospheric  model  adopted  (from  Machado  et  al.,  1980,  model  'FI'). 

Shown,  as  a  function  of  column  depth  N  ,  are  x  -  the  ionization 
fraction,  n  -  the  number  density,  T  -  the  temperature,  h  -  the 
height  above  the  photosphere  (t5qoo  *  D*  and  0  -  the  plasma 

resistivity,  assumed  classical  (s**"  Emslie  and  Machado,  1979). 

Figure  2.  Energy  deposition  rate  a  'solid!  and  depth-integrated  deposition 

rate  fQdZ  (dashed)  versus  N  .  Each  diagram  refers  to  a  different 

set  of  (v,  l  )  values  as  shown  around  the  border,  while  within 
o 

each  diagram  the  letters  a  through  d  refer  to  the  coronal  magnetic 

field  geometry  of  the  model  (see  Table  I  ) .  The  shaded  areas 

correspond  to  the  location  of  the  temperature  minimum  (cf.  Figure 

1);  Q  should  peak  in  this  depth  interval  for  the  model  to  be 

effective  at  heating  T  .  levels  of  the  atmosphere. 
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